Inverse Cascade Regime in Shell Models of 2-Dimensional Turbulence 
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We consider shell models that display an inverse energy cascade similar to 2-dimensional turbu- 
lence (together with a direct cascade of an enstrophy-like invariant). Previous attempts to construct 
such models ended negatively, stating that shell models give rise to a "quasi-equilibrium" situation 
with equipartition of the energy among the shells. We show analytically that the quasi-equilibrium 
state predicts its own disappearance upon changing the model parameters in favor of the establish- 
ment of an inverse cascade regime with K41 scaling. The latter regime is found where predicted, 
offering a useful model to study inverse cascades. 
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The inverse energy cascade in 2-dimensional Navier- 
Stokes turbulence is an important phenomenon with im- 
plications for geophysical flows In addition, it had 
been found that correlation functions and structure func- 
tions obey very closely Kolmogorov scaling (so-called 
K41), with only minute anomalous corrections, in con- 
tradistinction to 3-dimensional turbulence in which inter- 
mittency corrections to K41 scalingare sizable This 
difference is well documented j3[ but not yet under- 
stood. It is therefore tempting to construct simple mod- 
els of the phenomenon. Indeed, several attempts were 
made to construct shell models for this purpose, ^. 
So far these attempts ended negatively, failing to find a 
statistical steady state in which energy flows from smaller 
to larger scales together with having a Kolmogorov en- 
ergy spectrum. Rather, it was thought that whenever 
energy flew "backwards" , the statistical steady state set- 
tled close to thermodynamic equilibrium. In this letter 
we show that there actually exists a wide range of pa- 
rameter values for which shell models display the wanted 
behavior, thereby offering useful testing grounds for ideas 
on 2-dimensional turbulence. 

We discuss the issue in the framework of the Sabra 
shell model j^. Like all shell models |^ this represents 
a truncated Fourier representation of the Navier-Stokes 
equations. The Sabra model reads 
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where the dissipative term 7„ reads vk^ -\- ij,k~^f^ , with 
V and /i being the viscosity and drag coefficients respec- 
tively. Here u„ are complex numbers standing for the 
Fourier components of the velocity field belonging to shell 
n, associated with wavenumbers fc„. The latter are re- 
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stricted to the set /c„ — fcoA", with A being the spacing 
parameter, taken below to be 2. The forcing /„ is chosen 
here to act at intermediate values of n, n = rif, allow- 
ing in principle to study direct as well as inverse fluxes. 
The forcing is taken random with Gaussian time corre- 
lations as in the amplitude of the forcing is fixed 
below to l/A/2in all cases. The dissipative terms 7„ act 
both on the smallest and the largest scales with their re- 
spective (hyper)-viscosity and drag exponents a and (3; 
below we use a = (3 — 2. The dissipative terms become 
dominant at the viscous and drag scales rid and tll re- 
spectively. We will always have riL ^ nf <C rid- The 
coefficients a, h and c are adjustable parameters, with 
the constraint a -\-h -\- c = Q ensuring the conservation of 
energy in the dissipationless limit. Choosing a = 1 we 
explore the problem in terms of the single parameter 5, 
with -2 < 6 < §. 

It was shown before |l^ that for 6 < —1 there 
exist two positive definite invariants, the energy E and 
the "enstrophy" iJ, 
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which, in this case, are associated with an inverse and 
direct fluxes respectively However, the statistical 
steady state found in the regime —5/4 < h < —1 in 
[^, is close to thermodynamic equilibrium. This can 
be demonstrated via the properties of the structure func- 
tions, defined by 



S2{kn) = {\Un\^) , 

Ssikn) = Im{(w„_iu„M*_^i)} , 
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etc. Indeed, in Q these objects were found in the inertial 
range to be close to the exact solution in thermodynamic 
equilibrium which reads 
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Formula has two asymptotes: for small n in agree- 
ment with energy equipartition, and for large n with en- 
strophy equipartition. 
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Here ric ~ log(i3/^)/ log(a/c) is the cross over shell sep- 
arating the two asymptotic scaling forms of S'2. A and 
B are coefficients depending on the forcing and the dis- 
sipation. In particular, in this regime close to thermo- 
dynamic equilibrium, the cross over Uc moves to higher 
shells when the viscosity v is reduced. Unless otherwise 
stated, we choose parameters such that Uc^nf. 

Equation (0) implies zero fluxes. However in our sim- 
ulations we find in this regime a finite inverse flux of 
energy and a direct flux of enstrophy which do not go to 
zero when 7n ^ 0. The fact that the fluxes do not vanish 
also implies that ^3 is not exactly zero. One can write 
down the exact form of 6*3, which is correct always when 
there is a flux of energy or a flux of enstrophy : 



Ssikn) 



fc„ ^ , ul n Uf (energy flux), (11) 
-j , Uf <^ n nd (enst. flux) (12) 

A measure of the deviation of the statistics from Gaussian 
behavior is provided by the ratio 



which according to Eqs. (||)-(|l^) has the three separate 
regimes 
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These regimes are illustrated in Fig. |^. 

When R is small, it provides a measure of the magni- 
tude of the fluxes compared to their standard deviations. 
R is of order of unity at the dissipative boundaries, while 
it reaches its minimal value &t n — Uc- The former fol- 
lows from the fact that the dissipative boundaries are 
precisely where the second order dissipative terms bal- 
ance the third order transfer terms. In fact the ratio R 
cannot be larger than unity whenever scaling prevails. 
One sees this directly from the definitions ^ and (Q): 



Ssikn)/ \J S2(kri-l)S2{kn)S2{kn 



< 1 



(17) 



Since Uc moves to higher shells when the viscosity is re- 
duced, the value of R at the minimum decreases: we di- 
vide a decreasing S3 by an S'2 that remains constant over 
a larger range of n. We thus conclude that the quasi- 
equilibrium regime displays a alphabetical small param- 
eter when 1/^0. We will see that in the Kolmogorov 
regime there is only a numerical small parameter. 




Ssikn) / S2{k„y''^'^ as a function of 



*^ (dashed line) and u = 10"^^ 



FIG. 1: The ratio i?(fc„) 
n for b — —1.1 and 1/ = 10 
(solid line). One clearly sees the three regimes predicted by 
Eqs. (|l^, (p^, ([l6|), with the minimum occurring at n = n^. 
The minimum deepens and moves to higher shells when the 
viscosity is reduced. The vertical dotted lines designate the 
forcing scale. Note that at small values of n the statistics con- 
verges on much longer time scales than for large values of n; 
hence the remnant fluctuations around the scaling prediction. 



In ref. it was then discovered that there exists a 
transition for b crossing a critical value (b = —5/4 for 
A = 2) after which S'2 gains a new form in the direct en- 
strophy flux regime, close to the Kraichnan dimensional 
prediction ||l| 
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(up to small corrections). We note that this prediction 
can be inferred from Eq. ([l6| ) and the condition (|T^). 
Indeed R must be an increasing function of n towards its 
small scale boundary, which yields 
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or 6 > -5/4 for A = 2 and a = 1. Thus for b < -5/4 
Eq. ( |l6| ) can no longer be valid. While Ssffcn) docs not 
change, S2(kn) is replaced by the form (^8|) and conse- 
quently Eq. (|l^ is replaced by 
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In Fig. y we present this ratio as computed from numeri- 
cal simulations with the values of 6 = —1.5 and —1.6. We 
have used a total of 46 shells, with u = 10^'^^, /i = 10^^. 
The forcing was on shells 15 and 16. The three regimes 
are clearly seen, with the added important confirmation 
that this ratio is of the order of unity at the two dissipa- 
tive boundaries. 

Nevertheless previous work failed to find a similar phe- 
nomenon for the range of scales that supports the in- 
verse flux of energy. In that range the statistics remained 
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FIG. 2: Same as Fig. |l| but for b = -1.5 (dot-dashed line) 
and —1.6 (solid line). One clearly sees the three regimes pre- 
dicted by Eqs. (p^, (p^), (po|). Note that the constant regime 
is still effected by finite size effects from the viscous end. 
However the two curves are essentially the same. 

close to thermodynamic equilibrium, leading to the com- 
mon belief that shell models cannot be used to model 
2-dimensional turbulence. We explain next that the sta- 
tistical solution claimed for the regime b < —5/4, i. e. 
local thermodynamic equilibrium for the inverse flux of 
energy and direct enstrophy cascade, predicts its own 
destruction when b is reduced further beyond a critical 
value be that we can compute analytically. Indeed the 
set of Eqs. (||), (|l|), (|^) and the condition (0) fur- 
ther implies that R cannot be a decreasing function of n 
in the range ^ n <ti Uf, which implies 

(^)'^'>A4^6>-a(l + A-2/3) . (21) 

Accordingly, for b < be = — a(l -I- A~^/^) the quasi- 
equilibrium in the inverse energy flux regime can no 
longer be supported, and it changes into a true cascade 
regime with K41 scaling. For A = 2 this occurs at the 
critical value be ~ —1.63, where S'2(fcn) assumes the scal- 
ing form 

S2{kn) ^ ni, < n < . (22) 

Note that Eq. (|2^) implies the collapse Uc ul. In 
Fig. H we show the results of simulations at b = —1.9, 
with forcing at shells nf = 35, 36 and otherwise the same 
parameter values as in Fig. nl The agreement with K41 
scaling is apparent. We note that the scaling laws ( |2^ ) 
and (O) (which remains true in this regime) imply that 



R{kn) becomes constant as a function of kn- Thus we 
cannot display an alphabetical small parameter anymore. 
Nevertheless, the measurement of the constant value of R 
in the inverse cascade regime yields a number of the order 
of 0.02 or less. We thus have a numerical small param- 
eter, that is similar in magnitude to the corresponding 
value of R in 2-dimensional turbulence H . 




FIG. 3: Second (solid line), third (dashed line) and fourth 
(dot-dashed line) order structure functions in the inverse cas- 
cade regime for b = —1.9. The vertical dotted lines indicate 
the forcing range nf = 35,36. The dotted lines have K41 
slopes of -2/3, -1 and and -4/3 respectively. 

In summary, we exhibited a new regime of the statis- 
tical properties of shell models in which inverse energy 
cascade exists side by side with a direct enstrophy cas- 
cade. The statistical objects satisfy scaling laws in close 
correspondence with the Kraichnan dimensional predic- 
tions for 2-dimensional turbulence. Since this model is 
so much simpler than 2-dimensional Navier-Stokes equa- 
tions, it should provide useful grounds to understand the 
phenomenon theoretically. Such a discussion and a more 
detailed account of our numerical flndings will be pre- 
sented elsewhere p^ . 

Acknowledgments 

This work has been supported in part by the Euro- 
pean Commission under a TMR grant, the German Is- 
raeli Foundation, and the Naftali and Anna Backenroth- 
Bronicki Fund for Research in Chaos and Complexity. 
TG thanks the Israeli Council for Higher Education and 
the Feinberg postdoctoral Fellowships program at the 
WIS for financial support. 



[1] R. H. Kraichnan, Phys. Fluids 10, 1417 (1967). 
[2] P. Tabeling, Phys. Rep., in press. 



[3] L. M. Smith and V. Yakhot, Phys. Rev. Lett 71,352 
(1993). 



4 



[4] J. Parct and P. Tabeling, Phys. Fluids 10, 3126 (1998). 
[5] G. Boffctta, A. Colani and M. Vorgassola, Phys. Rev. E 

61, R29 (2000). 
[6] E. Aurell, G. Boffetta, A. Crisanti, P. Prick, G. Paladin, 

and A. Vulpiani, Phys. Rev. E 50, 4705 (1994). 
[7] P. D. Ditlevsen and I. A. Mogensen, Phys. Rev. E 53, 

4785-4793 (1996). 
[8] V. S. LVov, E. Podivilov, A. Pomyalov, I. Procaccia, and 

D. Vandembroucq, Phys. Rev. E 58, 1811 (1998). 
[9] T. Bohr, M. H. Jensen, G. Paladin, and A. Vulpiani, Dy- 



namical Systems Approach to Turbulence, (Cambridge, 

1998) and references therein. 
[10] E. B. Gledzor, Dokl. Akad. Nauk SSSR 209, 1046 (1973) 

[Sov. Phys. Dokl. 18, 216 (1973)]. 
[11] M. Yamada and K. Ohkitani, Prog. Theo. Phys. 79, 1265 

(1988). 

[12] T. Gilbert, V. S. L'vov, A. Pomyalov and I. Procaccia, 
in preparation. 



